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Abstract 

The Kramers- Wannier duality is shown to hold for all the even number spin correlation 
functions of the two dimensional square lattice Ising model in the sense that the high 
temperature (T > T c ) expressions for these correlation functions are transformed into 
the low temperature (T < T c ) expressions under this duality transformation. 

PACS: 05.70. Jk 

The two dimensional Ising model remains one of the very few models for which 
many properties can be calculated exactly' 1 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8,9 ]. Among these are the parti- 
tion function and in principle all the multi-spin correlation functions for the case of 
the square lattice (isotropic ) model. At least the algorithm for computing all the 
multispin correlation functions is known' 10 ' 11 ]. Quite a few of the two-spin correla- 
tion functions have been computed in terms of the complete elliptic integrals of the 
first and the second kind. It is also well known that the expression for the diagonal 
two-spin correlation functions for the anisotropic model have the same expression in 
terms of the complete elliptic integrals as that of the isotropic model. These corre- 
lation functions are known to satisfy a Painleve'VI nonlinear ODE with a proper 
definition of the variables' 13 ^ . All these things perhaps make it arguably the most 
widely studied model of ferromagnetism and antiferromagnetism 

Even though the model is so widely studied, new properties of this model are being 
discovered from time to time. It is the purpose of this communication to discuss one 
of such properties of the model. 

The very first sign that the model shows a phase transition from a disordered state 
to a ordered state was discovered by Kramers and Wannier ' 2 1 by considering what is 
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known as the duality property of the partition function of the model. We give here 
the Hamiltonian for the model which is 

H = -JJ2 WijVij+i + (!) 

n.n 

where o"jj = ±1, the products in the summation include only the nearest neigh- 
bours and i and j denote the positions of the spin on the lattice. The interaction 
strength J is positive for the ferromagnetic model and negative for the antiferromag- 
netic model. 

The partition function of the model can then be written as 

Z= E exp"^ (2) 

a=±l 

where/3 = (/c^T) -1 , and any multi-spin correlation function can be written as 

< Vi,jCTk,i ■••>=-= Yl a i,j a k,i ■ ■ ■ exp-^ (3) 

z a=±l 

The Kramers- Wannier duality can be understood by considering the high tem- 
perature and low temperature expansions of the partition function in the appropriate 
expansion variables. It is seen that each of the configurations of the high temperature 
expansion side has a counterpart in the low temperature expansion, and thus there 
is a one to one correspondence between the two phases of the theory. Kramers and 
Wannier' 2 ^ used this high temperature -low temperature duality to find the critical 
temperature for the model, the temperature at which the model undergoes a phase 
transition. According to their calculations, the critical temperature is given by 

(sinh2(3J) 2 = 1 (4) 



Later, in a celebrated paper, Onsager' 3 ! calculated the exact partition function 
by considering the transfer matrix for this model, which proved the correctness of 
Kramers- Wannier formula. The next quantity in line was the expctation value of a 
single spin, or in other words the spontaneous magnetization of the model, the form 
for which was given by Onsager and proved by YangW. Since then there were many 
numerical calculations of the multispin expectation values which finally resulted in 
writing down some of the two spin correlation functions of the model in terms of the 
complete elliptic integrals showing a beautiful structure of these functions when exact 
expressions are written for these correlation functions in terms of these integrals' 7 ' 8 ' 9 ! . 
It was also clear that all the two-spin correlation functions can indeed be written in 
terms of the complete elliptic integrals. The two-spin correlation functions were also 
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shown to satisfy certain quadratic difference equations that were derived by Perk' 12 ] . 
However it was not very clear how these quantities are related in the two phases of 
the model. 



In the meantime much work was done to understand the diagonal two-spin corre- 
lation functions and Jimbo and Miwa developed an equation that these correlation 
functions satisfy. General series solutions to this equation were dicussed in ref.' 14,15,18 '. 

For the exact expressions for some of the correlation functions of this type the 
reader is referred to ref.[7,17]. We can take any of these expressions as a function of 
the variable k = (sinh2/3J) 2 or 1/k — (sinh2/9J) 2 for T > T c ot T < T c respectively 
and consider the following non linear transformation on it. 

* I (5) 
K{k)->kK{k) (6) 
E(t) ^ (*) + (*■ -!)*(*) (7) 

K 

The last two transformation equations are just seen to be the correct transfor- 
mation rules for the complete elliptic integrals from their transformation properties. 
Under these transformations simultaneously, the high temperature (T > T c ) expres- 
sions are seen to transform into the low temperature (T < T c ) expressions and vice 
versa. This is true for all of the known diagonal two-spin correlation functions. How- 
ever there is nothing special about these correlation functions and this property must 
be true for all the spin correlation functions that can be written in terms of k and the 
complete elliptic integrals K (k) and E (k) and therefore this property must be true 
for all the even number spin correlation functions, which share this dependence on k 
and the complete elliptic integrals. It is essentially seen to be the Kramers- Wannier 
duality at work again. 

We now proceed to give a proof of the statement made above, from the Mon- 
troll, Potts and Ward's formulation' 5 ^ of even number spin correlation problem of the 
isotropic Ising model in terms of Pfaffians. In our notation, we will follow their paper 
as closely as possible. 

At first we consider the case of the rectangular lattice Ising model. Its partition 
function is given by 

Z= Y[e^p(pj 1 a i>j a i>j+1 + pj 2 cr itj a i+1;i ) (8) 

(t=±1 n.n 

where <7jj = ±1 is the state of the spin at the site of the latticed increases 
in one direction and j in the perpendicular direction. The product includes all the 
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nearest neighbours and the sum is over all the spins of the lattice. This, for the 
isotropic case J\ = J 2 reduces to expression(l). The transition temperature of this 
model is given by 

sinh 2/3 Ji sinh 2(3 J 2 = 1 (9) 

We introduce the variables 

z\ = tanh/3Ji z 2 = tanh/3J 2 (10) 
With the help of these variables the partition function can be written as 
Z = (cosh^Ji cosh(3J 2 ) N II + ^ij^ij+i) 

<r=±l n.n 

x (1 + z 2 a ijj a i+ltj ) (11) 

where N is the total number of spins on the lattice. The even number spin corre- 
lation functions can be formally written down as 

< ffo,oO" m ,^ P , 9 • • • >= Z^ 1 (cosh/3Ji cosh^Ja)^ x ^ cr 0:0 cr m:n a Ptq ■ ■ ■ 

a=±l 

X II U + WijVij+l) (! + WijVi+lj) ( 12 ) 

n.n 

where the product of spins on the l.h.s. contains an even number of spins. Following 
ref.[5] we can connect pairs of spins by (non overlapping) continuous line segments 
and we can introduce the factors of a^ k = 1 along those links and then arrange the 
resulting expression in such a way that it is related to the original partition function 
with a perturbed interaction strength along those links. One can then relate this 
expression to the Pfaffians by the usual correspondence. 

As discussed by the authors of ref. [5], this results in expressions which contain 
terms and integrals which can be related to the complete elliptic integrals of the first 
and the second kind in the case when we have the interaction strengths equal along 
both the vertical and horizontal directions, i.e, we have a square lattice. 

If we take the argument of these integrals to be the variable k in one phase of 
the system then the expressions for the other phase are obtained by simply replacing 
k by | everywhere in the expressions because there is only one expression for the 
correlation functions in terms of the Pfaffians, the one given by equation( 12). To make 
the argument clearer, we take one of the simplest cases, the case of < o"o,o°o,i > which 
is also one of the cases discussed in ref [5]. The generalization of the logic is obvious 
and therefore this case illustrates the properties of the general even number spin 
correlations quite simply. In the notation of ref. [5] this two-spin correlation function 
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is given by 

< o"o,o<7o,i 



>= (cosh/3Ji cosh/?J 2 ) 51 V 1 + z i 1(J o,oOo,i 

Z £T=±1 



x II ( 1 + ^l^ij'^ij+l) (! + z 2Vi,jCr i+1>:j 



(13) 



where the prime shows that the factors containing the product of spins which have 
occured previously in the expression have to be omitted in the later product. 

after going through some arguments this expression can be written as the product 
of the Pfaffian of 



z l - z x 



1 
-1 



and the Pfaffian of 



Zi 1 - Zx 





-i 



-[Zx -zx 



[0,1] 



RL 



[0,1] 



RL 







(14) 



(15) 



and zx, where [0, 1]^ is one of the elements of the matrix defined in the ref. [5]. The 
net result of all this is that 



< ^0,0^0,1 >= ± (zx - (l - zl) [0, 1] 



RL 



(16) 



This for the isotropic case whenzx = z 2 = z = tanh (3 J yields 



< ^0,0^0,1 >=±(z-(l- z 2 ) [0, 1] 



RL 



(17) 



However the z as well as the element [0, 1]^ can be written in terms of the variables 
k = {sinh2(3J) 2 or the low temperature expansion variable for which | = (sinh2/3J) 2 . 
As is clear the double integral expressions for the matrix elements involved can be 
written in term of complete elliptic integrals with the modulus k in the high tem- 
perature phase and modulus r in the low temperature phase. One therefore just has 
to look for the transformation properties of the complete elliptic integrals under the 
transformation 

and these are seen to be given by the equations (6) and (7)0. 



1 see for example ref. [16] 
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This therefore has provided us with a very important tool for relating the high tem- 
perature expressions to the low temperature expressions and vice versa. In general 
this means that we actually only need to know the T < T c expressions for any arbi- 
trary correlation function because the correlation functions involving an odd number 
of spins vanishes in the high temperature phase. Thus just the knowledge of T < T c 
expressions provide a complete description of the system. 

As a special case of these relations, in the Painleve'V I equation of Jimbo and 
Miwa, the solutions for integral values of n in the two phases are related to each 
other by the above mentioned transformations. But, of course, this remark is appli- 
cable only for the exact solutions and is not applicable to the series solutions of this 
equation term by term. 

We have thus learnt that the Kramers- Wannier duality has a greater validity than 
being applicable only to the partition function of the square lattice Ising model, it 
extends to all the even number spin-correlation functions of the model. 
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